Anomalous transport in Charney-Hasegawa-Mima flows 
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Transport properties of particles evolving in a system governed by the Charney-Hasegawa-Mima equation are investigated. 
Transport is found to be anomalous with a non linear evolution of the second moments with time. The origin of this anomaly 
is traced back to the presence of chaotic jets within the flow. All characteristic transport exponents have a similar value around 
fi = 1.75, which is also the one found for simple point vortex flows in the literature, indicating some kind of universality. 
Moreover the law ■y = fi + 1 linking the trapping time exponent within jets to the transport exponent is confirmed and 
an accumulation towards zero of the spectrum of finite time Lyapunov exponent is observed. The localization of a jet is 
performed, and its structure is analyzed. It is clearly shown that despite a regular coarse grained picture of the jet, motion 
within the jet appears as chaotic but chaos is bounded on successive small scales. 



I. INTRODUCTION 

Understanding transport in turbulent magnetized plasma is 
a task of fundamental importance. In these plasma, transport 
problems are often related to confinement, which is one of 
the last standing issues confronting the realization of mag- 
netically confined controlled fusion devices. As the literature 
evolves, there has been more and more evidence showing that 
the transport properties can be anomalous, in the sense that 
transport may not be correctly described by Gaussian kinet- 
ics, but by what one now calls "strange kinetics" |JJ] . In these 
regards, transport phenomena in turbulent magnetized plasma 
is part of the ever growing number of physical systems dis- 
playing anomalous properties I2l Ij, 0, y[ . As of today the full 
understanding of these phenomena are far from being com- 
plete and in many regards a full blown theory able to capture 
and describe correctly these "strange kinetics" has not yet sur- 
faced. There seems nevertheless to be a common agreement 
to link these phenomena to Levy-type processes and their 
generalizations. Moreover the use of fractional derivatives 
in Fokker-Plank-Kolmogorov type equations captures quali- 
tatively some of the transport properties and is thus a good 
step towards a proper description of anomalous transport |6]. 

The link between the Hamiltonian dynamics and the kinet- 
ics at origin of anomalous transport properties are relatively 
well understood when dealing with low dimensional systems 
such as a time periodic flow which belong to the class of 
1 — 1/2 degree of freedom Hamiltonians. The dynamics in 
these systems is not ergodic: a well-defined stochastic sea, 
with chaotic dynamics, filled with various islands of quasi- 
periodic dynamics compose the phase space. The anomalous 
properties and their multi-fractal nature are then linked to the 
existence of islands within the stochastic sea and the phe- 
nomenon of stickiness observed around them However 
when dealing with more complex systems the loss of time- 
periodicity complicates the picture. For instance in geophys- 
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ical flows or two-dimensional plasma turbulence, the islands 
which were static and well localized in phase space, are re- 
placed by "coherent structures", which have a life of their 
own. Hence, tackling the origin of anomalous transports from 
the chaotic dynamics of individual tracers becomes more sub- 
tle. Recently, the existence of a hidden order for the trac- 
ers which exhibits their possibility to travel in each other's 
vicinity for relatively large times was exhibited |9]. This or- 
der is related to the presence within the system of chaotic jets 
fTolfTlll . These chaotic jets can be understood as moving clus- 
ters of particles within a specific domain for which the motion 
appears as almost regular from a coarse grained perspective. 
Typically, the chaotic motion of the tracers is confined within 
the characteristic scale of a given jet, within which nearby 
tracers are trapped. From another point of view, looking for 
chaotic jets can be understood as a particular case of measure- 
ments of space-time complexity | ll|- 

The purpose of this paper is to study transport properties 
and to look for chaotic jets in a model of two-dimensional 
turbulence which applies either in the plasma context where 
it is known as the Hasegawa-Mima equation 11311 or in the 
geophysical one where one speaks about the Charney equa- 
tion (see for instance lfl4lo . In this setting Annibaldi et al 
have already shown strong evidence of anomalous transport 
of passive particles I13l Hal , hence the search for chaotic jet 
is expected to give some clues on the origin of anomalous 
transport in these systems, and for instance to identify the 
structures responsible for such transport. Indeed in |9], it 
was clearly shown that trapping within chaotic jets resulted 
in anomalous transport. Moreover it appeared that jets were 
localized around the coherent structures of systems, namely 
the vortex cores and that the structure of the jets itself was a 
hierarchy of jets within jets, reminiscent of the multi-fractal 
nature of transport observed in 1 — 1/2 degree of freedom 
Hamiltonians systems |8]. Hence, the goal of this paper may 
be seen as two-fold. First we want to understand the origin 
of anomalous transport in a model of two-dimensional plasma 
turbulence. And second, since we are using chaotic jets to 
track this origin, we are at the same time testing the existence 
of these jets in a more complex setting than the system of point 
vortices used in \9]. Once the presence of jets is confirmed we 



2 



may be able to speculate that the different anomalous transport 
behavior portrayed in the nonexhaustive following references 
iHEoEDjllSi] may all find their origin in the presence 
of long lived chaotic jets in the considered systems. 

The paper is organized as follows. In Sec. |H]the basic def- 
initions are introduced. A brief introduction of the Charney- 
Hasegawa-Mima equation is given. Then the dynamics of test 
particles is given and the Lagrangian approach of transport in 
this setting is discussed. In Sec. [TI]]the dynamical evolution of 
the field and of passive test particles, as well as transport prop- 
erties are computed numerically. First the numerical setting 
is discussed. Then three different regimes corresponding to 
different choices of the parameters for the field evolution are 
considered cases. In Sec. UVI we track the origin of anomalous 
transport in the three considered. First we recall the definition 
of a chaotic jet and present the numerical method used to de- 
tect these jets. We then present the statistical results related 
to trapping time within jets and analyze the origin of anoma- 
lous transport by localizing the jets and by describing their 
structure. Finally we conclude in Sec. [V] 

II. BASIC DEFINITIONS 

A. The Charney-Hasegawa-Mima equation 

The Charney-Hasegawa-Mima equation can be written in 
the following form, 

d t Q =0 (1) 

fl = <P-\A$ + gx, (2) 

where [•, •] corresponds to the Poisson operator, ft is a gener- 
alized vorticity given by Eq.0, $ is a potential and A and g 
are parameters. 

Typically Eqs. Q and (0 can either describe the evo- 
lution of an anisotropic plasma, and are then referred to as 
the Hasegawa-Mima equation, or the evolution of geostrophic 
flows in which context they are known as the Charney equa- 
tion. This formal identity has an advantage as the results ob- 
tained in this paper for transport properties should apply in 
either context. 

It is however important to be able to put the results back 
in a physical context. With this in mind we provide in the 
next section a short derivation of Eqs Q and Q in both the 
anisotropic plasma configuration and the geostrophic approx- 
imation. 



is adiabatic, n e (x,y,t) = no(x,y)e e ^^ x ' v ' t ^ T,! , where n e is 
the electron density, no is the equilibrium plasma density, and 
T e the electron temperature. In the anisotropic plasma let the 
cyclotron frequency beuj c = eB/rrii, and let the hybrid sound 
speed be c s — y/T e /rrii, where is the ion's mass. 

We now consider the motion of cold ions (the ion temper- 
ature is assumed zero, Ti = 0) on characteristic time and 
length scales much larger respectively than l/u c and the De- 
bye length \d- In this situation the plasma is quasi-neutral 
fii ~ n e , and a linear combination of ion continuity equation 
and momentum equations gives (see 12111 for details) 

„ D { e<j>\ „, 1 D . 

-V.v_l = — — +v ± .Vlnn = — ■ -prrfw+Wc) , 

Dt \T e J oj + oj c Dt 

(3) 

where vi is the speed in the plane perpendicular to the mag- 
netic field and uj = V x vj_ is the vorticity of the 2D flow. The 
slow motion of the ions necessarily implies that vj_ ~ \r D = 
-Bq~ 2 Bo x V</> and, using the notation (p = ecj)/T e , subtraction 
of the right hand side terms in Q readily gives 

— (A^-^/i|)-v D .Vlnn = 0, (4) 

where 1^ = c s /lj c is the hybrid Larmor radius. 

When the plasma is homogeneous (no=const.), the tur- 
bulent dynamics of @ is characterized by the absence of 
waves and the formation of large vortices. Moreover, in the 
II — * +oo limit, Eq. @ becomes formally equivalent to the 
2D momentum Euler equation (jp being the stream function in 
the Euler case). 

In plasma devices, density is larger in the core than at the 
boundaries. It this setting we can often write no(x,y) ~ 
noe _a; / L ™, where L n is a characteristic density gradient 
length. Using this last approximation we finally obtain the 
Hasegawa-Mima equation: 

£(Z*-WS)-4S-0, (5) 

where g = u> c /L n . 

The inhomogeneous character of the equilibrium density 
profile in Eq.(|5} implies the existence of so called drift-waves 
in the flow. In particular, the waves deform and interact with 
the self-organized vortex-like or multipole-like structures, and 
play a key role in their interactions. For instance the collision 
of two dipoles can lead to two monopoles plus one dipole plus 
some radiation l22ll . 



B. Wave-vortex paradigm equation for two-dimensional flows. 

1. The Hasegawa-Mima equation 

Let us start by considering a magnetized anisotropic 
plasma, e.g a plasma with a uniform magnetic field along a 
direction, B = Bz. We shall also assume that electron re- 
sponse to the turbulent fluctuations of the electric potential cf) 



2. The Charney Equation 

Let us now explain the physical nature of the intrinsic 
length in atmospheric motions. In the context of shallow water 
approximation where the density p is uniform, the modeling 
of the evolution of the atmospheric winds leads to a dynamic 
equation similar to or (|5jl- Indeed, the atmosphere is char- 
acterized by small Ekman and Rossby numbers 1 14], so that 
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the motion of thin incompressible fluid rotating layer (thin 
with respect to the characteristic scale L± of the horizontal 
motion) is such that friction forces are very weak compared to 
the inertial Coriolis force, second and that pressure forces are 
balanced by gravity in the vertical direction z (hydrostatic ap- 
proximation). As a consequence the horizontal pressure gra- 
dient is independent of the vertical component z. It is also 
proportional to the layer width h(x, y, t) (because of the free 
boundary condition p(x, y, h) = po = const.). It is therefore 
reasonable to assume that the horizontal velocity field vj_ is 
also z-independent: v = vi(i, y, t) + w(x, y, z, t)z. 

If we assume that the rigid earth surface is locally flat, i.e 
w(z = 0) = 0, integration of the incompressibility condition 
leads to w — — zV.v^and, in particular, dh/dt = — /iV.vj.. 
Moreover, since the forces acting on the fluid are the Coriolis 
and pressure forces, pfv± x z and — pS/p, it is straightfor- 
ward to deduce from the momentum equation the exact iden- 
tity d(u + f)/dt = (lu + /)V.vj_. It follows that 



C. Advection equation 

For an incompressible fluid, the evolution of a passive par- 
ticle is given by the advection equation 

z = v(z,t) (8) 

where z[t) represent the position of the tracer at time t in the 
complex plane, and v(z, t) is the velocity field. An important 
feature of this evolution of passive particles is that the evolu- 
tion equation given by Eq. can be rewritten as: 

a$ t .0$ 

z = ~l-^= , z = i— , (9) 

oz oz 

where the potential $ acts as a time dependent Hamiltonian, 
and~denotes the complex conjugate. This Hamiltonian struc- 
ture is fundamental as it imposes some constraints on the dy- 
namics of passive tracers, which should be taken into account 
when carrying a numerical simulation. 



V.Vj. 
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Dt 



h 
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TfDi 



(w + /) , (6) 



where we put h = h + Hq, Hq being the mean width of the 
layer and h the fluctuations around the mean (h <C Hq by hy- 
pothesis). The Coriolis term fz is the local component of the 
planetary vorticity uj p i anet in the vertical direction, the north- 
ward local component f n y being negligible. Indeed, the Cori- 
olis force in the horizontal direction is {—2pio v i anet x v)j_ = 
P v ± x / z + P wz x fnY ~ pv± x fz because, in the shallow 
water approximation, w/v± = 0(Ho/L±) <C 1 and far from 
equator latitudes, /„ ~ /. 

It is clear from equations and (jfji that in the special 
limits where the plasma is homogeneous (hq =const.) and 
the Coriolis parameter / = fo is taken constant, both phys- 
ical systems are formally equivalent. We just have to iden- 
tify e</>/T e with h/Ho and choose as length and time units 
(l L ,l/u} c ) or (Ijt, 1/ /o)> where l R = \/gH a /f a is the so- 
called Rossby length. But the similarity holds even for the 
case of a layer with a slow dynamics (u> <C /) by taking into 
account the slow variation of the Coriolis parameter in the 
northward y-direction: / = fo + Py with Sf /fo <C 1 (/3-plane 
approximation). In fact, since the Rossby number is small, we 
can also make the approximation ~ vjj = gf~ 1 z x Vh. 
And by subtracting the right hand side (r.h.s) terms of (j§Jl, we 
obtain the so-called Charney equation 



DV 



(7) 



It is clear that @ and Q are formally equivalent. However, 
in the atmospheric case, waves are not generated by a density 
inhomogeneity but by the earth rotation. As mentioned earlier, 
both equations Q and Q can be written in the compact form 
given by Q where Eq. actually writes CI = <p — l^Aip — 
x/L n (the transformations tp h, g (3, and Ir lead 
to the Charney equation). 



D. Numerical settings 

1. Charney-Hasegawa-Mima 

Simulations of the Charney-Hasegawa-Mima equation are 
performed for different initial conditions and choices of pa- 
rameters. The choices are made to be consistent with the lit- 
erature, namely the conditions chosen in j3l fl^l . The simula- 
tions are performed within a square box of size L = 20 and 
periodic boundary conditions using a pseudo-spectral code. 
In order to compute the evolution of passive tracers accu- 
rately we settled for a somewhat low resolution mesh of 128 2 . 
Fourier transforms are computed using a fast Fourier algo- 
rithm. For the time evolution, we chose a 4th-order Runge- 
Kunta integration scheme with typical time step St = 0.05. 

In order to avoid numerical instability, as well as a trivial 
asymptotic behavior Eq. Q could not be kept as it is and a 
dissipation term D as well as a forcing term F were added: 



d t Q + \Q, $1 = D + F . 



(10) 



One may also argue that the dissipation term may be relevant 
to describe some physical phenomena. 

For the dissipation we used a hyper-viscous term D = 
^(Vf2) 4 . For the forcing we used a term, whose Fourier trans- 
form is: 

1 



F(k x ,k y )=F Q 



(11) 



where for any (k x ,k y ), (p(k x ,k y ) is a random phase uni- 
formly distributed on the circle; k runs over an intervall cen- 
tered on ko with a range 5k. One can set also a time de- 
pendency by regularly updating the realization of the random 
phases. 

The initial condition is given by the following choice for 
the SI' = $ - AA$ field: 



ft' = A ' 



1 
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L 



2im 



L 
(12) 



y + <Pi,j 



where ipij is a random phase uniformly distributed on the cir- 
cle. 



2. Passive tracers 

It is important to take special care of the way the dynam- 
ics of the tracers are computed to characterize the possible 
anomalous properties of transport. If they exist, such anoma- 
lies should find their origin in the existence of "memory ef- 
fects", namely long time correlations, since the accessible 
range of velocities is finite. In this perspective, any source 
of randomness leading to memory loss due to the numerical 
scheme may then induce a spurious effective diffusive behav- 
ior. Moreover the Hamiltonian nature of the tracers dynamics 
imposes necessarily the choice of a simplectic integrator. 

We thus chose the sixth-order implicit Gauss-Legendre 
simplectic scheme to compute the trajectories [23], as this in- 
tegration scheme was successfully used in systems of point 
vortices |3 l25ll . However in order to avoid a possi- 
ble source of noise, we had to compute the velocities of par- 
ticles "exactly", meaning that we performed an exact back- 
Fourier transform of the modes describing the evolution of 
the field. This constraint is numerically expensive, and ex- 
plains our choice of a relatively low resolution of 128 2 for 
the evolution of the field. In this setting the evolution of pas- 
sive tracers may be understood as describing the advection of 
particles in a flow field generated by 128 2 modes interacting 
through Eq. Q. 



E. Field settings 
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Table I: Density of Energy and Enstrophy for the three considered 

cases 





Energy 


Enstrophy 


Smooth Field 


0.64 ±0.01 


0.8 ±0.05 


Forced Field 


2.62 ± 0.005 


10.81 ±0.03 


Anisotropic Field 


0.28 ± 0.005 


0.18 ±0.01 
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Figure 1 : Visualization of the field — A<3? for the choice of parameters 
A = 1, g = 0.1, v = 7 10" 6 , L = 20, N = 128 2 with no forcing. 
The field is "smooth" and appears as being isotropic. A few vortices 
are present. The gray coloring scheme scales from -4.5 (black) to 4.5 
(white) 



As to the choice of parameters and initial conditions in 
Eq. we considered three different cases, with quite differ- 
ent values of the parameters. As mentioned earlier we chose 
as a starting base similar initial conditions, type of forcing, 
and size of the system as in 1 3 ] . For the three considered cases, 
for given initial conditions, we let the system evolve until it 
becomes sufficiently stationary for the time span considered 
during simulation, Tfi na i = 10 4 . Stationarity is considered as 
being reached by monitoring the evolution of the energy and 
enstrophy of the system (see Table|I}. When forcing is present 
(case 2 and 3), the random phases of the forcing are updated 
all at once every St = 2. The values assigned to Ao are re- 
spectively for the three cases: 1.1, 0.001, 0.5. For the first 
case, as there was no forcing, we tried a decent amplitude, for 
the second case with strong forcing we let the forcing drive the 
system and for the last case, we tried something intermediate. 
The motivation for these choice was to get three different sta- 
tionary fields, as the focus of this paper is actually more the 
transport of test particles rather than the actual dynamics of 
the field. 

In order to visualize the field we chose to use levels of the 
function — A$. The three different considered cases are rep- 
resented in figuresm2]and|3] 



1. Smooth Field 

To obtain the "smooth" field depicted in Fig. ^ we carried 
out simulations with no forcing (i*b = 0) and a low dissi- 
pation. The parameters for this run were A = 1, g = 0.1, 
v = 7 10~ 6 . Due to this low dissipation the energy may be 
considered as constant for the length of the simulation. We 
notice that a few distinct vortices are present. During the evo- 
lution a merger between two vortices occurred. There is also 
an average drift in the y-direction. 



2. Forced Field 

To obtain the "forced" field depicted in Fig. [2] we carried 
out simulations with a strong forcing and dissipation, the pa- 
rameters for this run were A = 4, g = 0.1, v = 5 10~ 5 , 
Fq = 4, ko = 6 ± 2, the phases for the random forcing were 
updated every At = 2 time units. The value of ko corre- 
sponds to physical scales of Sx ~ I 21 « 1. With this choice 
of parameters the system consists of two perturbed vortices. 
An average drift in the y-direction is also noticeable. 
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Figure 2: Visualization of the field — A$ for the choice of parameters 
A = 4, g = 0.1, v = 5 10~ 5 , F = 4, k = 6 ± 2 , At = 2. Two 
big perturbed vortices are present. The gray coloring scheme scales 
from -1.4 (black) to 1.4 (white) 



Figure 3: Visualization of the field — A$ for the choice of parameters 
A = 0.125, g = 2, v = 7.5 10~ 6 , F = 1.5, k = 12 ± 2, At = 2. 
Elongated structures in the y-direction are present, marking a strong 
anisotropy. The gray coloring scheme scales from -1.6 (black) to 1.6 
(white) 



3. Anisotropic Field 

To obtain the "anisotropic " field depicted in Fig. |5] we 
carried out simulations with some forcing and a high value 
for g. The parameters for this run were A = 0.125, g = 2, 
v = 7.5 1CT 6 , F = 1.5, k = 12 ± 2. The value of k 
corresponds to physical scales of Sx ~ |p ~ 0.5. The phases 
for the random forcing are updated every At = 2 time units. 
In this settings elongated structures as well as a strong drift in 
the y-direction is present. 



III. TRANSPORT PROPERTIES 

A. Definitions 

Unfortunately the deterministic description of the motion 
of a passive particle in a chaotic region is impossible. Lo- 
cal instabilities produce exponential divergence of trajecto- 
ries. Thus even an idealized numerical experiment is non- 
deterministic, as round-off errors are creeping slowly but 
steadily from the smallest to the observable scale. The long- 
time behavior of tracer trajectories is then necessarily studied 
by using a probabilistic approach. In the absence of long-term 
correlations, the kinetic description, which uses the Fokker- 
Plank-Kolmogorov equation, leads to Gaussian statistics. Yet 
if a phenomenon with associated long time correlations oc- 
curs, profound changes in the kinetics can be induced. These 
memory effects sometimes result in the modification of the 
diffusion coefficient in the FPK equati on 12a. 12711 . But often 
their influence is more profound 1 7, 8, 9, 28, 29t l30ll . and leads 
to non-Gaussian statistics, and for instance to a non-diffusive 



behavior of the particle displacement variance: 

(( S -( S )) 2 )~^, (13) 

where (• • •) stands for ensemble averaging. Within this proba- 
bilistic approach, the main observables in order to characterize 
transport properties will be moments of the distributions: 

M q (t) = (\s(t) - ( S (tW) , (14) 

where q denotes the moment order. The finiteness of velocity 
and of time in our simulations implies that all moments are 
finite and that a power law behavior is expected 

M q ~ D q t^ q) . (15) 

with, generally, ^ q/2 as would be expected for normal 
diffusion. The nonlinear dependence of /i(<?) is a signature of 
the multifractality of the tr ansp ort, while its linear dependence 
reflects a fractal situation IDOm 1 L 13 211 . In the fractal situation 
all of the moments can be described by a single self-similar 
exponent v 

= v ■ q , (16) 
whereas the case when /j,(q) is nonlinear 

plq) = ^ cons t (17) 

q 

transport is multifractal. This distinction is important since in 
the weak case the PDF evolves in a self-similar way: 

p(M) = e = t-"(*-<*)) (is) 
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while a non-constant v{q) in precludes such self- 

similarity (see the discussions in |71 @] f° r details about the 
non self-similar behavior). 

In order to characterize transport, we focus on the arc length 
s(t) of the path traveled by an individual tracer up to a time t, 

Si(t) = f Vi (t')dt' , (19) 
Jo 

where Vi(t') is the absolute speed of the particle i at time t'. 
This choice is motivated by the fact that in order to consider 
mixing properties from the dynamical principles it is impor- 
tant to consider the trajectories within the phase space. How- 
ever, for the system of passive particles the phase space is for- 
mally "identical" to the physical space where particles evolve 
(see Eq|9jl- Another important feature of the arc length is that 
it is independent of the coordinate system and as such we can 
expect to infer intrinsic properties of the dynamics. More- 
over, expression d!9l > implies that S{ (t) reflects the history up 
to time t of the speed Uj. Since the velocities are bounded 
the observations described further on of anomalous transport 
behavior are directly linked to strong memory effects. 



B. Particle Transport 

We followed the evolution of 512 passive particles as de- 
fined by Eq.ljSJ, and we computed Sj(i)'s up to a time r = 
10 4 . In order to keep a constant numerical accuracy the in- 
crements Asj(nT) = s,((n + 1)T) — Si(nT) were recorded 
for the successive diagnostic times nT. This last feature has 
also the advantage of providing better statistics when comput- 
ing moments. Indeed since we have chosen stationary regimes 
for the field, we can assume that transport properties are inde- 
pendent of the initial condition of the field and thus use time- 
invariance to increase statistics. This time invariance becomes 
however not useful for times close to t — 10 4 , and statistics 
in this region is thus not accurate anymore, also we are not in- 
terested in short time behavior. We have therefore considered 
times 10 2 < t < 5.10 3 in order to compute the fi(q). 

In all three considered cases the transport is found to be 
anomalous and superdiffusive, with a characteristic second or- 
der exponent /i(2) w 1.8 for all cases. The time behavior of 
the moments and characteristic exponents for all considered 
cases are plotted in figures @J [5] and|S] and a summary is pro- 
vided in table|IT| where we included for comparison the results 
obtained in Ref. |9J] for a flow governed by point vortices. 
The similarity observed for the exponents in these quite dif- 
ferent settings of the parameters and regimes of the Charney- 
Hasegawa-Mima equation as well as the one observed in point 
vortices, points to a universal behavior for the transport of pas- 
sive tracers in these two-dimensional flows, which in some 
sense confirms the validity of the quite general estimation of 
/i fts 1.5 proposed in 1 8] as a value for an universal exponent. 

However, conversely to the point vortices case, one can no- 
tice from Figs. E]|5] and [6] that the behavior of characteristic 
exponent versus moment for the first two cases may corre- 
spond more to the multifractal type although the nonlinearity 




q 



Figure 4: Top: Moments of distribution of the arc length M q (r) = 
(|s(t) — (s(t)) \ q ) versus time of tracers evolving in the smooth field 
for q = 1/2, 1, 3/2, • ■ ■ 8. The behavior M„(r) ~ r" (g) is con- 
firmed. Bottom: Characteristic exponent versus moment order, q 
vs fj,(q), with fi(2) — 1.87. o and + symbols refer respectively to 
computations made using 512 and 256 particles. The solid line cor- 
responds to fi(q) — (J,(l)q expected for self-similar behavior. Trans- 
port is super-diffusive and multifractal. 



is quite weak, and error bars are quite large. For the third case, 
transport is more or less simply fractal. This behavior corre- 
sponds to simple fractal transport in the anisotropic system, 
which implies an almost self-similar behavior of the distribu- 
tion function. In fact as may be illustrated in figure|5j we may 
expect that with more particles and larger times the multifrac- 
tal behavior may be more clear, a feature which can also be 
the case for the third case. When computing characteristic ex- 
ponents one also has to recall the possibility of log-periodic 
oscillations [33]. This phenomenon may also be responsible 
for the uncertainties on the measured values of the exponents. 

All in all these results imply that at least for intermediate 
times, transport is anomalous and single fractal, which means 
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Figure 5: Top: Moments of distribution of the arc length M q (r) = 
(|s(t) — {sir)) \ q ) versus time of tracers evolving in the forced field 
for q = 1/2, 1, 3/2, • • • 8. The behavior M q (r) ~ t m(9) is con- 
firmed. Bottom: Characteristic exponent versus moment order, q vs 
fi(q), with fi{2) = 1.68. One can notice a change of slope after 
log 10 (ro) = 3.3. The o and + symbols refer respectively to compu- 
tations made using 512 and 256 particles for r > to, while the * sign 
corresponds to r < To. The solid line corresponds to = [i(l)q 
expected for self-similar behavior. Transport is super-diffusive. 



Figure 6: Top: Moments of distribution of the arc length M q (r) = 
(\s(t) — (s(t))| 9 ) versus time of tracers evolving in the anisotropic 
field for q = 1/2, 1, 3/2, ■ ■ ■ 8. The behavior M q (r) ~ r^ {q) is 
confirmed. Bottom: Characteristic exponent versus moment order, q 
vs fj.(q), with /i(2) = 1.85. o and + symbols refer respectively to 
computations made using 512 and 256 particles. The solid line cor- 
responds to fi(q) — (J,(l)q expected for self-similar behavior. Trans- 
port is super-diffusive and single fractal. 



Table H: Characteristic second moment exponent for the three dif- 
ferent cases studied. Exponents obtained in flows governed by point 
vortices in Ref. 1 9] are given for comparison. 



Point vortices 


4 vortices fi(2) » 1.82 
16 vortices n(2) » 1.77 


Charney-Hasegawa-Mima 


Smooth Field ft(2) ta 1.81 
Forced Field ^(2) » 1.73 
Anisotropic Field ft(2) m 1.85 



that transport properties in these systems should be correctly 
described by a fractional Fokker-Plank-Kolmogorov equation 
of the type [6]: 



dPP(s,t) d a P{s,t) 



= V- 



d\i 



(20) 



with fj, sa 2/3/a (see for instance |6, 8]). For larger time trans- 
port may remain single fractal or develop a multifractal behav- 
ior in which case a model of transport properties using Ea.i20i 
may be only approximate. Note that using Ea. J20> to model 
transport properties implies that a kinetic limit has been per- 
formed on particle statistics. In this situation we are dealing 
with Levy type processes, hence moments higher than two are 
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Figure 7: Tracking of e coarse-grained regular jet. 



likely not defined. This particularity is however linked to the 
kinetic limit, and since we are dealing numerically with a fi- 
nite number of particles during a finite time, and also since 
velocities are bounded, moments of particle arc-length distri- 
bution are defined, see 1 3411 for a discussion on the matter. 



IV. JETS 
A. Definitions 

Tracking the origin of anomalous transport can be fairly 
well understood when one is able to draw a phase portrait 
using a Poincare map and for instance to measure Poincare 
return time to a given region of phase space. The conclusion 
of this type of analysis will almost certainly lead to the fact 
that the phenomenon of stickiness on the boundaries of the is- 
lands generates strong "memory effects", as a result of which 
transport becomes anomalous. However, when dealing with a 
more complex system, for which the drawing of a phase por- 
trait is not achievable, one has to rely on other techniques. 

In a two dimensional phase-space, the phenomenon of 
stickiness corresponds often to passive particles remaining 
for large times in the neighborhood of an island of regular 
motion. A consequence of this behavior is that sticky zones 
are regions where particles are trapped, and therefore are re- 
gions where particles remain in each others neighborhood for 
large times. It becomes then natural when dealing with more 
complex systems for which no phase portrait can easily be 
drawn, to look for places where passive particles remain in 
each other's vicinity for large times. One possibility to cap- 
ture this feature of the dynamics is to look for its signature 
by measuring finite time Lyapunov exponents (FTLE), and by 
isolating within the space of initial conditions, regions of van- 
ishing FTLE's (see for instance 1 35 ]). This type of approach 
has however its shortcoming, namely sticky regions are not 
necessarily smooth from the microscopic point of view, mean- 
ing they can be regions of strong chaos that are somehow re- 
stricted within an arbitrary small scale, which may be prob- 
lematic when dealing with FTLE's. The other possibility is 
to look directly for chaotic jets 0|. These chaotic jets can 
be understood as moving clusters of particles within a specific 
domain for which the motion appears as almost regular from a 
coarse grained perspective. From another point of view, look- 
ing for chaotic jets can be understood as a particular case of 
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Figure 8: Tail of the distribution of trapping times r for the smooth 
field (Fig.EJ. A power-law decay is observed with typical exponent 
p(r) ~ r" 7 with 7 fa 2.8 ± 0.1. 

measurements of space-time complexity 1 12). 

In order to look for jets we proceed as described in the il- 
lustration presented in Fig. Qwhich provides an easy and in- 
tuitive description of the mechanism used to detect jets. To 
summarize, we consider a reference trajectory r(t) within 
the phase space. We then associate to this trajectory a cor- 
responding "coarse grained" equivalent, i.e the union of the 
balls U.B(r(t), e) of radius e whose center is the position r(t). 
Given an e-coarse grained trajectory, we analyze the behavior 
of real trajectories starting from within the ball at a given time, 
and measure the time r and length s, before the trajectory es- 
capes from the coarse grained trajectory. We then analyze the 
resulting distributions. This approach has already been used 
with success when studying numerically the advection of pas- 
sive tracers in flows governed by point vortices |9]. The main 
difficulty in using this diagnostic follows from the fact that 
data acquisition is not sampled linearly in time nor space, a 
point which leads in the present case to some difficulties. 



B. Statistical results 

In the setting of the evolution of passive tracers within the 
three considered Charney-Hasegawa-Mima flows, we settled 
for the following values e = 10 _1 and 5 = 10~ 3 . These 
values are to be compared with the e = 10~ 2 and S = 10~ 6 
considered for point vortex systems in |9). This reduction of 
accessible scales stems from the fact that for the length of our 
simulations (t = 10 4 ), only a few trajectories are escaping 
from the jets when we try to use the point vortex values. This 
leaves us with not enough data to gather realistic statistics. It 
is important to note that results should not depend much on the 
value of e (see the discussion in |9]), as long as its value does 
not cross a characteristic scale in the system. For instance 
in the point vortex systems, cores surrounding vortices had 
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Figure 9: Tail of the distribution of trapping times r for the forced 
field (Fig.|2j- A power-law decay is observed with typical exponent 
p(r) ~ t -7 with 7 « 2.8 ± 0.1. 
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Figure 10: Tail of the distribution of trapping times r for the 
anisotropic field (Fig.[3}. A power-law decay is observed with typical 
exponent p(r) ~ t~ 7 with 7 » 3.5 ± 0.5. 



a radius of order ~ 0.1, or here in the forced systems, the 
characteristic wavelength corresponds to scales of order ~ 1. 
Moreover in order to infer some possible fractal structure of 
the jet, we are constrained to consider the largest possible e 
once 5 has been fixed. Here the choice made for 5 was the 
smallest we could make in order to gather enough statistics, 
although these may not be sufficient for the anisotropic case. 
We also tried to keep at least two orders of magnitude between 
e and <5, to infer eventual fractal properties. 

With the present choice of parameters we were able for 
gather about 15000 events for each of the three cases. We 
followed 256 reference trajectories for a time r = 10 4 dur- 
ing which the behavior of two nearby tracers was checked. 
We call a reference trajectory the trajectory of a given passive 



tracer which evolves freely, we then put randomly at a distance 
6 of this tracer a second (or more) tracer which was called a 
ghost in 1 9]. We then let these tracer evolve until the distance 
e is reached. Then the ghost is removed and time interval At 
and travelled length As are recorded. A new ghost is assigned 
to the reference trajectory and so on. We of course expect that 
the portion we compute of reference trajectory (r = 10 4 ) is 
sufficiently ergodic in the accessible phase space. 

We were then able to obtain the trapping time distributions 
/?(t) described in Figs.[SJ[5]andE3 The characteristic expo- 
nent p(r) ~ t 1 for trapping times observed in two out of 
the three systems is typically 7 1=3 2.8. For these cases we 
therefore obtain good agreement with the 

7«M(2) + 1, (21) 

relation which links the transport exponent p, to the character- 
istic trapping time exponent. The law Hit can be linked to 
the fractional transport equation i20\ . One can link the 7 ex- 
ponent to the (5 exponent as follows (3 = 7— 1 , and p = 2/3 /a 
(see 1 8, 9] and references therein), a being linked to the spa- 
tial fractal properties, the equation d2 11 1 is valid when a ~ 2. 
One can also notice in Fig. s. |SJ [5] that the law p(r) ~ t 1 
applies typically for times up to t ~ 1000 or a little more. 
This behavior can be expected as we are following a finite 
number of reference trajectories (256) during a finite time 
(Tmax = 10 4 ). In order to potentially obtain a broader range 
of applicability one, would have to increase the value of r max 
to at least 10 5 . This is unfortunately beyond our computing 
resources. Note also, that due to the inverse cascade there is 
an accumulation of energy at large scales, such that the driving 
field may not always be considered to have reached a station- 
ary state if we decided to increase simulation length. In order 
get better statistics, one also has to be careful in not adding 
to many particles instead of increasing time. Indeed the limits 
N — > 00 and r — > 00 are most likely not commuting in these 
anomalous regimes where rare events related to memory ef- 
fects are important. 

In fact the law ( 12 li also applies in the anisotropic field for 
small trapping times, but o. verth whome range it is clear in 
Fig. ^3 m at no scaling law really emerge. When comparing 
with the two other cases, we may find two possible reason 
for this failure. First, we can notice that in this configuration 
of the field the energy and enstrophy are quite low, imply- 
ing weak nonlinear effects. It may thus take more time to 
fill the tail of the trapping time distribution. Second, another 
peculiarity of the anisotropic case is the the fact pointed out 
in l3rj l37il : in the anisotropic case, the conservation of the 
generalized vorticity ft = $ — AA$ + gx . by the evolu- 
tion of the Hasegawa-Mima equation Q implies necessarily 
that for strong values of g the motion along the x— direction is 
bounded. Hence in this situation the motion of passive tracers 
is quasi one dimensional. In this setting it is likely that the 
hidden fractal properties in the system may differ than for a 
two-dimensional systems, especially regarding the derivation 
of the laws linking (3, 7 and a. In these regards, the fact that 
the law Mil is valid for small times up to a time r c may just be 
a simple consequence of the fact that for r < t c , the passive 
tracers have not yet reached the boundaries imposed by g and 
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Figure 11: Unnormalized Distribution of Lyapunov exponents od, 
see Eq. 1231 . An accumulation towards zero is observed. The zoom 
near small regions reveals an actual decrease for an < 0.02, such 
behavior is expected as data is bounded by finite speed and finite 
time. 



a spatially two dimensional behavior is still accurate. At this 
point it is important to recall that we have periodic boundary 
conditions. The constraint imposed by the conservation of fl 
is thus not as important for small values of g corresponding to 
the first two cases. 

The rarity of events with smaller values of 5 indicates that 
the trajectories of tracers are relatively regular when one con- 
siders small scales, which implies the possibility of vanishing 
Lyapunov exponents. Hence, following the definitions given 
in 11: 

° L = i ln i' (22) 

°d = -r- In | , (23) 
As o 

we compute the distribution of Lyapunov exponents <jq ob- 
tained from the smooth case data in Fig.^2> where one can 
directly see an accumulation towards zero of the distribution. 
The zoom in Fig. II llshows an actual decrease for cd < 0.02. 
This behavior is expected as we only have finite speeds in the 
system and simulations are carried for a finite time. For in- 
stance do < 0.02 corresponds to an average trapping time 
t a = 1100. This accumulation of exponents toward zero is 
quite important as it is a strong evidence of memory effects. 
The system does not display universal hyperbolic properties, 
and thus transport modeling using this as a hypothesis are 
probably doomed. This also shows the relevance of the no- 
tion of weak chaos and weak mixing properties, and therefore 
the necessity of considering simple models such as the billiard 
considered in 1 38]. Indeed these provide good insights on how 
the presence of quasi zero Lyapunov exponents does not nec- 
essary mean the system is regular, but that complexity is not 
growing as fast as expected 11 1 211 . 
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Figure 12: Localization of a long lived jet in the "Forced Field" (cf 
Fig.|2j- The jet is bouncing back and forth between the two perturbed 
vortices 



C. Localization and structure of jets 

Due to the shape of the finite size Lyapunov exponent dis- 
tributions observed, it is also possible to track and localize 
jets which are responsible for the anomalous behavior. In- 
deed, these exponents are monotonically decreasing. Hence 
we can define a threshold beyond which the jet can be con- 
sidered "regular". We can track the reference trajectory after- 
wards in order to localize regions responsible for anomalous 
behavior. In Fig.^]we show the localization of a jet for which 
the trapping time of nearby tracers is found to be w 1000, in 
the case of the forced field (cf Fig.|2j. The jet is bouncing back 
and forth between the two perturbed vortices. The procedure 
to detect a jet is quite simple, we chose a reference trajectory 
and start to compute trapping times of ghosts. When we notice 
that a ghost has not escaped for a while (we choose r > 200 
for this jet), we start recording the positions of the reference 
tracer. We also reset the ghosts and add many more in order 
to track the "structure" of the jet, and record their positions as 
well. When one of the ghosts leave the jet, we measure the 
actual trapping time. If it appears long enough (here we chose 
t > 1000), we stop, else we reset everything and wait until 
the reference trajectories detect a new potential jet. The ad- 
vantage of this method is that we have no a priori information 
on the jet location. For instance the jet depicted in Fig. ^] 
was not anticipated, as we expected more a trapping within a 
vortex or in its neighborhood. We may notice however that 
by using this method of detection we miss the first portion of 
the jet. Note also that by localizing jets, we are able to detect 
regions responsible for anomalous behavior of transport prop- 
erties. It therefore may therefore give a good clue on where 
and how to act on a system, if one wants to reduce or eliminate 
this anomaly. 

Once a jet is located, one can infer the behavior of the test 
tracers with respect to the reference trajectory. In order to get 
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Figure 13: Relative position of test particles for the long lived jet de- 
picted in Fig. ll2l ftop). Zoom (bottom), test particles reach distances 
< 10 -4 to the reference trajectory. 



Fig- El 256 test particles were initialized on a circle around 
a reference particle whose trajectory was trapped in the long 
lived jet localized in Fig. El ar >d the relative positions of the 
test particles were plotted during the life of the jet. One no- 
tices a "star" like shape of the figure, this phenomenon is a 
consequence of multiple stretching and squeezing of the cir- 
cle of particles, which one can observe while looking dynam- 
ically at the evolution of the circle. One also can notice that 
ghosts which initially are localized at a distance S = 10~ 3 
from the reference tracer, can find themselves much closer to 
the test particles close to 10~ 4 in Fig. El This implies that 
we do not expect that the value of the trapping time exponent 
7 is strongly dependent on the initial choice of S. Note also 
that the behavior portrayed in Fig. El seems to indicate that 
chaos is present at these scales, thus contrary to our physi- 
cal intuition the relative motion within the jet is not regular. 
However despite this fact, the chaotic behavior of trajectories 



is trapped within a given small scale for a long time, which in 
the end gives rise to a regular non-diffusive trajectory from a 
coarse grained perspective. This trapping phenomena is also 
illustrated by the apparently "quantified" maximum radii seen 
on the star shape displayed in Fig. El We may thus speculate 
that these radii correspond to actual successive barriers which 
are restricting chaos within a small scale. Note also that the 
chaotic behavior within the jet may have induced numerical 
computations of strong values for finite time Lyapunov expo- 
nents. In this situation relying on these exponents would have 
meant that jets which are responsible for anomalous transport 
may go unnoticed. 



V. CONCLUSION 

In this paper we have investigated the dynamical and statis- 
tical properties of passive particle advection in different con- 
figurations of Charney-Hasegawa-Mima flows. The goal of 
the work was to consider transport properties of these sys- 
tems while putting in perspective the results obtained for point 
vortex flows. In this sense it was a step further in provid- 
ing qualitative insights on general transport properties of two- 
dimensional flows. The transport properties of all the con- 
sidered cases are found to be anomalous with characteristic 
exponent fj, ~ 1.75 — 1.8. These values are also quantitatively 
comparable to the results obtained for point vortex flows. 

In order to analyze the origin of the anomalous transport 
properties, passive tracer motion is analyzed by measuring the 
mutual relative evolution of two nearby tracers, i.e by look- 
ing for chaotic jets |0]. The jets can be understood as mov- 
ing clusters of particles within a specific domain where the 
motion is almost regular from a coarse grained perspective, 
inducing memory effects and long time-correlations. The dis- 
tribution of trapping times in the jets shows a power-law tail 
whose characteristic exponent is in very good agreement with 
the law 7 = [i + 1 linking the transport exponent /i to the 
trapping time exponent 7. This agreement is a good signature 
that the origin of anomalous transport in these system is inti- 
mately related to the existence of jets in the sense described 
previously. The localization of jets in the system can be done, 
and it is shown that jets are not necessarily located around a 
coherent structure as was the case for point vortex flows, but 
that they can manifest themselves by a trajectory bouncing 
back and forth between 2 structures. Moreover when analyz- 
ing the "structure" of the jet, it is shown that the trajectory 
of individual tracers are likely chaotic within the jet, but that 
this chaotic behavior is for long times restricted within a given 
small scale, giving rise to the regular non-diffusive structure. 
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